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Abstract 

Recent experiments reveal both passive subdiffusion of various nanoparticles and anomalous active transport of such 
particles by molecular motors in the molecularly crowded environment of living biological cells. Passive and active 
microrheology reveals that the origin of this anomalous dynamics is due to the viscoelasticity of the intracellular fluid. How 
do molecular motors perform in such a highly viscous, dissipative environment? Can we explain the observed co-existence 
of the anomalous transport of relatively large particles of 100 to 500 nm in size by kinesin motors with the normal transport 
of smaller particles by the same molecular motors? What is the efficiency of molecular motors in the anomalous transport 
regime? Here we answer these seemingly conflicting questions and consistently explain experimental findings in a 
generalization of the well-known continuous diffusion model for molecular motors with two conformational states in which 
viscoelastic effects are included. 
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Introduction 

After the publication of Albert Einstein's tlieory of Brownian 
motion in 1905 [1], Jean Perrin reported the first systematic 
microscopic studies of individual diffusing particles in 1908 [2]. 
Today, modern single particle tracking techniques routinely reveal 
insight into the stochastic motion of submicron tracers in aqueous 
solutions at unprecedented resolution, thus allowing one to drrectiy 
observe the transition from initial ballistic to diffusive Brownian 
motion [3] and to resolve the effects of hydrodynamic backflow 
[4] . Measuring the passive and driven motion of microprobes has 
become a standard means to characterize soft matter [5]. 
Particular attention is currendy paid to the relaxation and 
diffusion dynamics in dense colloidal systems [6,7] and inside 
hving biological cells [8,9]. 

The intracellular fluid (cytosol) of biological cells is a superdense 
[10] aqueous solution containing biomacromolecules such as 
proteins and RNA at volume fractions of up to 40%, a state often 
referred to as macromolecular crowding [11,12]. Indeed, the state 
of crowding in the cytosol effects severe changes of the diffusion 
behavior of submicron particles [8,13,14]. Thus, anomalous 
diffusion of the form <,[dt{t)]^} = i[t(t+t')-t{l'fy^tl' with 
0 < /? < 1 is observed in living cells for the passive motion of single 
biopolymers, endogenous granules, viruses, and artificial tracer 
particles [15-21]. Compared to normal Brownian motion with 
P=l, these particles therefore subdiffuse [22]. Such anomalous 
dynamics presents a challenge to the development of controlled 



uptake of drugs and nanoparticles and their intracellular delivery 
by molecular motors for therapeutic processes [23]. 

What happens to the active motion of particles in the cytosol 
which are driven by molecular motors [24—32], see Figure 1? WiU 
the overall dynamics of the coupled motor-cargo system be 
affected by the superdense state of the cytosol, and how? In fact, 
anomalously fast diffusion was observed in living cells for the active 
motion of viruses, microbeads, and endosomes [16,17,33]. When 
the motors normally transport their cargo with some mean velocity 
in a given direction, ballistic superdiffusion with /J = 2 is measured 
[16]. However, various subbaUistic power exponents l<p<2 are 
also found [17,33]. Here we come up with a mesoscopic physical 
approach for molectilar motors in the cytosol of living cells and 
show that our predictions are in agreement with the experimen- 
tally observed behavior. The simplest explanation for the 
experimental fact that transport by motors along intracellular 
microtubuli yields subbaUistic superdiffusion (1 <j8<2), is provid- 
ed by the sublinear scaling {3x(t)y = <[x(t + r') — x(t')ycct''''" of 
the mean position x along a single microtubule, with an effective 
scaling exponent 0 . 5 < 0£eff < 1 . Then, assuming quenched disorder 
in microtubuli orientations without a preferred direction yields 
/? = 2o£j.ff which will range in the subbaUistic superdiffusive regime. 
Can normal and anomalous transport with various p co-exist and 
be mediated by the same motors in the same cell? Which role is 
played by the size of the cargo, and what determines the precise 
behavior of such active transport and its efficiency? These are the 
questions we answer in this article. 
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Figure 1. Molecular motors moving along a microtubule in the 
crowded cytoplasm. A large cargo is subject to viscoelastic drag, 
effecting dramatic changes in the transport dynamics. 
doi:1 0.1 371 /journal.pone.0091 700.g001 

A well-established physical approach to anomalous transport 
phenomena is based on the intrinsic viscoelasticity [5,34—37] of 
complex fluids such as the cytosol. Depending on its size and 
speed, a nanoparticle may feel different effective viscosities in the 
macroscopic limit of normal diffusion which may be enhanced by 
a factor of hundred to several thousands with respect to the 
viscosity of pure water [38^0]. Transient regimes of anomalous 
diffusion become biologically important when the corresponding 
spatial length scale is comparable to the linear size of cells, 
typically several micrometers. The viscoelastic nature of the 
cytosol and other dense solutions has been verified in several 
experiments [19,20,41,42]. 

Results 

Ratchet Model of Molecular Motors 

We study the interplay of the viscoelastic environment of the 
cytosol with the action of a molecular motor and its cargo. A well- 
established model of Brownian motors of the kinesin family is 
based on the continuous diffusion of a Brownian particle in a 
potential landscape, which randomly fluctuates in time between 
two realizations, V\(x)^V2(x), depending on the internal state of 
the motor, that is undergoing active conformational fluctuations 
[25-29,31,32]. These conformational fluctuations are caused by 
the binding of negatively charged ATP molecules to the motor 
(state 1) and the reactions of ATP hydrolysis and dissociation of 
products (ADP and phosphate group) making up state 2 within a 
minimalist modeling framework. The potentials V\ 2(x) describe 
the free energy profiles leading the motor molecule in the 
corresponding conformational states along a microtubule. Since 
microtubules are periodic dipolar structures with period L « 8 nm, 
the potentials reflect this periodicity, Vx jix-'r L) = Vi 2{x). More- 
over, Fi(a: + L/2)= V2{x) within our two-state motor model, such 
that two potential switches occur during one cyclic turnover of the 
motor enzyme (power-stroke or "hand-over-hand" mechanism) 
and advance the motor by one period L. The motor direction is 
determined by the polarity of the microtubules, reflected in 
the space inversion asymmetry of the potentials. We use the 
harmonic mixing model [32] with Vi(x) = — Uq cos {2nx / L) — 
U\ sin(47tx/L), see the upper inset in Figure 2. This potential has 
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Figure 2. Normal transport for large cargo particles, large 
potential amplitude and small turnover rate, in the absence a 
of constant loading force, /o = 0. Single motor transport (full line) is 
almost perfectly locked to the potential fluctuations (broken red line 
depicting a renewal process counting the number of potential 
fluctuations in units of L/2) occurring with mean turnover frequency 
i'tuni = 112 Hz, in a potential (top inset) with amplitudes Uq = 0.25 eV 
((7o=l in dimensionless units) and f/i =0.162 eV (AK = 0.7 eV), for 
L = 8 nm. A particle with an effective radius 300 nm (like a magnetic 
endosome [33]) experiences asymptotically for r»T,„iu5 ~22.4 sec an 
effective viscous friction enhanced by a factor of jjctr = lO'' with respect 
to water. The bottom inset shows that on the relevant transient time 
scale the free particle subdiffuses with anomalous diffusion coefficient 
Z), *368 nm^/s" Initially, diffusion is normal. The time-average over a 
single trajectory, <<5.y^(0>7-= 737 Jq'^ ' [x(t + t')-x(t')]^dt' , is shown for 
r = 0.1 sec and compared with the theoretical subdiffusive ensemble- 
averaged result (red line). See iVIethods. 
doi:1 0.1 371/journal.pone.0091 700.g002 



a metastable state within each period. Thermal fluctuations play a 
positive role, allowing the motor to avoid getting trapped in such 
metastable states on its search for the potential minimum after 
each conformational change due to ATP binding and ATP 
hydrolysis. 

Within the power-stroke idealization, the maximal mean 
velocity of the motor due to the fluctuations between the potentials 
V\;i(x) becomes v = Lvturn, where Vtum is the typical motor 
turnover rate. The latter is composed of the conformational 



transition rates Vi_2 according to v. 



-1 . 

turn ' 



+ .As the simplest 



approach based on the widely applied flashing ratchet model 
[26,29,31,32,43], we assume that these transition rates do not 
depend on the transport coordinate x. Such an ideal motor would 
consume one ATP molecule (with an energy amount of 0.5 to 0.62 
eV or 20 to 25 ksT at room temperature, ks being the Boltzmann 
constant) while transferring a cargo over the distance L. The 
corresponding energy E\j^(t) invested into the temporal increase of 
the potential energy of the motor during repeating turnover cycles 
of the ^catalytic wheel' [44] can be calculated as a sum of potential 
energy jumps A V(x(ti)) occurring at random instants of time ti 
marking cyclic conformational transitions Ki F2 -> Ki -> . . . [45] . 
Furthermore, apart from delivering a cargo over a distance &x{t), 
the motor can perform useful work, W^!x(t)=fQ5x(i), against 
some constant force opposing its direction of motion. Ensemble 
averaging over many trajectory realizations, we obtain the 
thermodynamic efficiency i^th = < f^^use(0)/^-'^m(0)- the long 
time limit this is a time-independent quantity in the normal 
transport regime, where both and <(fFuse(0) ^^re 

proportional to time. If ^0=0, the thermodynamic efficiency is 
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zero, and all the input energy will be dissipated as heat. In order to 
characterize the energetic performance of a molecular motor in 
such a situation, we introduce the delivery efficiency D defined as 
the ratio of deliver)' distance d to the product of delivery time / 
and the average number of turnover cycles (Afturn) (number of 
ATP molecules consumed). D thus has the meaning of a mean 
delivery velocity per input energy amount (in dimensionless units). 
The goal is to deliver cargo over a certain distance as quickly as 
possible using the smallest amount of energy. Thus, ideally 
</ = L<iVtum> and <iVtura> = Vturn'- Therefore, Adeal = 
d/{t<^Ntumy) = L^Vi^m/d, which hnearly increases with the 
turnover frequency at some fixed d. However, we expect the 
delivery efficiency to deviate from this idealization. 

Both the motor and its cargo are subjected to friction and 
random thermal forces from the environment. For normal viscous 
Stokes friction, the frictional drag force is /fr(/)= —t]Qx{r), where 
f/o is the friction coefficient. It is proportional to the medium's 
viscosity f and the particle size. For a sphere of radius a, rjQ = 6naC. 
If we assume that the linker between the motor and its cargo is 
rigid, we can model their coupled motion as that of a point particle 
moving under an effective frictional force. In this way, one 
accounts for the cargo size simply by adjusting the effective 
friction. The dynamics for such a simplified motor is then defined 
by the Langevin equation 



nox(t)=- 



dV(x,t) 

8x 



-fo + (o(t). 



(1) 



where <^o(0 is a Gaussian random thermal force with zero mean, 
which is completely characterized by its autocorrelation 
function. The fluctuation-dissipation relation •(io(0^o('')) = 
2kBTr]gd(t — t') ensures that the description is compliant with 
the laws of thermodynamics, so that no directional motion can 
emerge when V(x,t) is time-independent and fixed to either Vi{x) 
or V2(x). The prototype motor model (1) has been investigated in 
great detail within a Markovian setting, using various model 
potentials and a different number of motor substates [25- 
27,29,32,43]. 

Viscoelastic Environment 

The above ratchet model is appropriate to describe the motor 
dynamics at dilute solvent conditions in vitro. Our focus here is 
different, as we want to study the motor action in living cells, 
where the following experimental facts have been established: 
Even in the absence of cargo the effective friction coefficient for 
the motor is enhanced by a factor of 100 to 1000 in the cytosol 
compared with the one in water [28,38]. This phenomenon is due 
to the superdense state of the cytosol, crowded with various 
biopolymers. Concurrently, the diffusing nanoparticles themselves 
experience a medium with an effectively enhanced viscosity, that 
can exponentially depend on the particle size [39,40]. Moreover, 
numerous experiments reveal [20,33,36,37] that the complex 
shear modulus G'iw) of the cytosol displays a power law scaling 
G*((b)oc (/(»)", a ranging between 0.2 and 0.9 for frequencies in 
the range from inverse milliseconds to several hundred inverse 
seconds [5,9,20,33,37]. This reflects the viscoelastic nature of the 
cytosol, which needs to be taken into actxjunt for molecular motors 
[17,33,34,46-51]. To explicitly consider viscoelastic effects for 
cargo particles is even more pressing given the experimental results 
revealing viscoelasticity-induced subdiffusion of free, passively 
moving submicron particles [17,20,21,41,42], whose size is 
comparable to typical, larger cargo such as vesicles. This implies 
that: (i) the dynamics must be described by a frequency-dependent 
friction corresponding to a viscoelastic memory for the friction 



term with a power-law kernel [34] ((f) ocf such that (ii) the 
diffusion of free particles becomes anomalously slow (^dx^(t)yccf 
with 0 < a < 1 on the corresponding time scales from milliseconds 

to minutes. This is a mesoscopic, transient effect. However, it 
becomes very important for transport processes in living cells as 
they occur on exactly the physiologically relevant time and length 
scales. Consequentiy, the Langevin equation (1) for the motor's 
velocity must be extended to the generalized Kubo-Langevin form 



. . , dV(x,t) . „ . . 



n„,(t-t')x{t')dt'^(ijt). 



(2) 



where r]„(t)ccl^{t), and imit) represents colored thermal 
Gaussian noise with autocorrelation function <^m(0^m(^)/' = 
kgTrji^dt — fl), as demanded by thermodynamics (absence of 
directed transport in a static periodic potential K,(x)) and the 
Kubo second fluctuation-dissipation theorem [52]. On physical 
grounds, a memory cutoff always exists for rj„{t) such that on a 
sufficiently long intermediate time scale the kernel has the scaling 
property V]^{f)azt^^ . This ensures that the effective friction 
coefficient f/o + fJeff with )7efr= lo^ ^m(')d! is finite, but strongly 
enhanced over tJq, in compliance with previous studies of 
molecular motors [28]. There exists also a short time cutoff, 
corresponding to the largest vibrational frequency of the medium 
contributing to the friction. Given the upper and lower cutoff, one 
can approximate the memory kernel t]„,(t) by a finite sum of 
exponentials [34,53] and the generalized Langevin equation (2) 
can be derived from a corresponding multi-dimensional Markov- 
ian Maxwell-Langevin model of viscoelastic dynamics [34], see 
Methods for details. On the basis of this theoretically and 
experimentally well-founded approach the results herein were 
obtained from numerical analysis. 

Perfectly Normal Transport 

The first major surprise is that even carr\'ing large cargo 
particles like magnetic endosomes with radius about 300 nm our 
model motor can operate by an almost perfect power stroke 
mechanism in the normal transport regime, as demonstrated in 
Figure 2. For this, the binding potential amplitude should be 
sufficiently large, AF = 0.70 eV, and the turnover rate sufficientiy 
small, Vturn —100 Hz, both reasonable values for this motor-cargo 
system. This almost perfect dynamics occurs in spite of the fact 
that the free cargo alone subdiffuses, <i5x^(?))-2Z)a/^/r(l +£>:), on 

a transient time scale between Tmax 

/'Jeff' (see Methods) 

and Tmax with a subdiffusion coefficient as small as 

v=22.4 sec, and 



Z)„«3.68x lO-K-mViO-^ for j;, 



eff- 



:10l 



a = 0.5, an experimentally relevant value we consider in this 
work. This normal transport regime is possible as viscoelastic 
subdiffusion is ergodic [8,34,53]. Moreover, for a sufficiently large 
potential barrier AV separating spatial periods of the potential the 
particle has enough time Trei to relax and settie down sufficientiy 
close to the potential minimum. Thus, it can be advanced by a 
half-period at each switch of the potential, despite the power-law 
character of the relaxation dynamics, if only the time scale 
separation condition Trei « Tturn « Tesc is satisfied. This condition is 
easy to fulfill for realistically small turnover frequencies v and 
realistic AV-15 to 30 fc^TVoom since the mean escape time grows 
exponentially fast with the barrier height, Tosc^^ sxpiAV/kgT). 
The occurrence of such normal transport is consistent with most 
observations on molecular motors. Thus we accomplished our first 
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goal, to explain the active normal transport of particles, that 
otherwise subdiffuse when they are not attached to a motor. 
However, the model motor considered here is actually assumed to 
be a little too strong. It has a stalling force of about 1 5 pN, almost 
twice the typical value, see below. Let us therefore consider a 
somewhat weaker binding potential and faster enzyme turnover 
rates to see the origin of anomalous transport for such a large 
cargo. 

Anomalous Transport 

Indeed, when we decrease the potential amplitude by 2/3, and 
consider an increasing turnover frequency in the presence of an 
opposing (external drag) force _/o, which further reduces the 
effective barrier height, an anomalous transport regime 
<(5x(0>ocf"="' with 0<aeff<l is enforced, see Figure 3,a. The 
thermodynamic efficiency in this strongly anomalous transport 
regime is rather small, see Figure 3,b. Moreover, it becomes an 
cdgebraiccJly decaying function of time [50,51], Rthif)'xi/t^~'^ ■ 
This happens because Wi,^c(r)cC'(5x(t)yccf''' scales sublinearly 
with time, while the input energy' is consumed at a constant rate. 
This means that asymptotically most of the input energy is used to 
overcome the dissipative influence of the environment character- 
ized by the massively enhanced elfective viscosity. Concurrentiy, 
the useful work performed against the force fo always remains 
finite and the stalling force is also about the same as for normal 
diffusion ratchets. This finding agrees with topical experimental 
values of 7 to 8 pN for the staUing force of kinesin motors [28]. 
The slower the motor turnover is, the larger the effective transport 
exponent (Zcff becomes, along with a higher motor efficiency. The 
efficiency displays almost a parabolic dependence on the force fo, 
Rth'xfoi^ —fo/fo) in this regime. A similar dependence was 
derived for fluctuating tilt viscoelastic ratchets [50,5 1] . The stalling 
force is not only roughly proportional to the barrier height, but 
also depends on the frequency of the potential switches. 

Smaller Cargo Size 

Depending on the cargo size and the binding potential 
amplitude the transport can become more normal and thermo- 
dynamicaUy highly efficient even for a large turnover frequency, as 
Figure 4 illustrates for Vtum = 1.12 kHz. This is about the maximal 
turnover frequency which can be expected for molecular motors. 
For this, the potential amplitude should be sufficiently large and 
the cargo smaller in size. Here, we reduced fj^.^ to 10^ with the 
same Tmax = 22.4 sec. Hence, assuming that the effective viscosity 
of the medium remains unchanged, '7c(-'7cff'^mM' becomes reduced 
by a factor often, which corresponds to a cargo with one tenth of 
the size, that is, of some 30 nm radius. However, if we were taking 
into account its dependence on the particle size [38-40], this value 
of ^jfj should in fact be attributed to a somewhat larger particle. 
The subdifTusion coefficient is enhanced accordingly, 
Z);, «3.68T0^" m^/s"-^. Furthermore, to show that the fractional 
friction strength rj^ and the subcliffusion coefficient are 
characteristic for the transport properties rather than rj^^ and 
Tmax separately, we also considered the case with fjef[ = W^ and 
Tmax = 2240 sec yielding the same D^, see data with Uo = \.0 in 
Figure 4. For the largest potential amplitude in Figure 4 the 
thermodynamic efficiency is appreciably high, up to 45% for the 
studied case. This is very surprising: The transport efficiency in the 
anomalous regime can be temporally almost as high as the 
maximal efficiency of kinesin motors in the normal regime (about 
50%). For this potential amplitude, however, our motor is stronger 
than a typical kinesin motor. It has a staUing force of about 1 5 pN, 
see Figure 4. The effective exponent a^s is about 0.87 at the 
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Figure 3. Anomalous transport of large cargo particles at lower 
potential amplitude, larger turnover rates, and in the presence 
of loading force fo. (a) Effective anomalous transport exponent a^t 
and (b) thermodynamic efficiency i?th while working against a constant 
force fo near the end point of the simulations (0.224 sec or 10^ in 
dimensionless units). The thermodynamic efficiency decays over time as 
^tii(0°c l/'' """- The analysis considers the same particles as in Figure 2, 
but here the potential height is reduced by factor of 2/3. Ensemble 
averaging is performed over 10' particles and random realizations of 
potential flashes. The inset in (a) shows the dependence of izeff on the 
mean enzyme turnover frequency for fo = 0. 
doi:1 0.1 371/journal.pone.0091 700.g003 



maximum corresponding to « 10 pN. However, the transport 
is anomalous and the efficiency decays algebraically as 
i«th(Ooc?-° ". For increasing loading force fo the anomalous 
diffusion exponent a^{[ becomes smaller and the thermodynamic 
efficiency drops faster as function of time. This means that the 
optimal value of the force fo corresponding to the maximum of 7?th 
slowly shifts towards smaller values, as if the motor became 
gradually ever more tired, and more quickly exhausted for a 
higher load. Upon reduction of the barrier height by a factor of 
3 /4 the transport is still close to normal for = 0- Moreover, the 
thermodynamic efficiency can still be temporally rather high at 
optimal load. However, (Zeff now drops faster with fo. A stronger 
reduction of the potential amplitude, by one half, leads immedi- 
ately to the emergence of a low efficiency, strongly anomalous 
transport regime, even for /o = 0 (Figure 4). Still, a strong 
reduction of the turnover frequency down to Vtum = 100 Hz wiU 
recover the normal transport regime for a sufficientiy small _/() in all 
cases considered. The transport of even smaller particles is clearly 
normal for realistic parameters. 

Delivery Performance 

For vanishing loading force fo = 0 thermodynamic efficiency is 
zero, even in the normal transport regime. We proposed above 
that the performance of molecular motors such as kinesin should 
be characterized by the energetic efficiency D of the cargo delivery 
over a certain distance d. We calculated the above-defined 
delivery efficiency D in Figure 5 as function of the turnover 
frequency Vmrn for several different deliver)- distances d, potential 
amplitudes Uo, and fractional friction coefficients i;^. Remarkably, 
for a small cargo (smaller fj^ff in Figure 5) the calculated 
delivery efficiency follows the ideal power-stroke dependence, 
D = L?'Vxa^ld in the entire range of realistic turnover frequencies. 
Even for a relatively large cargo. Figure 3, but for much smaller 
turnover frequencies the transport is close to this ideal normal 
regime. It is expected to be normal already for Vtum = 10 Hz. Such 
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Figure 4. Dependence of (a) the effective transport exponent 
otefl and (b) tlie thermodynamic efficiency Rti, on the load fo, for 
three different potential amplitudes and for turnover frequen- 
cy Vtiirn = 112 kHz. Ensemble averaging is done over 10^^ particles and 
random realizations of potential flashes, »;cff = 10', Tniax = 22.4 sec, or 
f/(,iT=10'' and Tninx = 2240 sec (with the same see Methods). 
Matching of the results for two sets with the same amplitude (Jo = 1.0 
indicates that i/^, is the characteristic quantity, rather than ))(,,-(■ and r^ax 
separately. Efficiency is calculated at the end point of simulations 10^ 
corresponding to ? = 0.224 sec. 
doi:1 0.1 371 /journal.pone.0091 700.g004 

a small frequency is, however, not accessible to numerical analysis. 
In the anomalous regime, with increasing turnover frequency the 
delivery efficiency reaches a maximum at v^ax and then decreases. 
The anomalous transport becomes less efficient for v > v^ax and 
Vmax shifts to smaller values with increasing delivery distance d. 
We speculate that in the cell economy this effect can become 
relevant for the optimization of the metabolic budget. 

Discussion 

We proposed a simple basic model which reconciles experi- 
mental observations of both normal and anomalous transport by 
highly processive molecular motors in biological cells. Our model 
presents an immediate generalization of a well-known two-state 
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Figure 5. Delivery efficiency D as function of the turnover 
frequency Vtun,. In the ideal power-stroke regime, Adcai = i^i'turn/<^- 
For small Vtum our results agree with this simple dependence (broken 
lines). 

doi:1 0.1 371 /journal.pone.0091 700.g005 



model of normally diffusing molecular motors by accounting for 
the viscoelastic properties of the intracellular ffuid. It not only 
explains how molecular motors may still operate by a power-stroke 
like mechanism while carrying a large cargo which subdiffuses 
when left alone, but also why and how an anomalous transport 
regime emerges for even larger cargo. It is crucial for this 
explanation that viscoelastic subdiffusion and anomalous transport 
possess finite moments of residence times in any finite spatial 
domain. Thus, there exist time scales for sliding down towards the 
potential minimum within one period of the flashing ratchet 
potential, for the escape to another potential well, and for the 
mean turnover time of the potential flashes. When the time scales 
are well separated, the transport is normal. However, when the 
sliding time scale starts to interfere with the turnover time, or the 
interweU potential barrier is lowered so that backsteps can occur, 
the transport becomes anomalous. These qualitative basic features 
are expected to survive in more complex models of molecular 
motors operating in viscoelastic environments. 

Specifically, we showed that transport by molecular motors 
becomes anomalous for large cargo particles with large fractional 
friction coefficient rj^ when the enzyme turnover is fast, and the 
binding potential amplitude AK is not sufficiently large. Larger 
potential amplitude A K for a fixed spatial period leads to faster 
relaxation of the motor particle to a new potential minimum after 
each potential switch and thus to more normal transport behavior. 
The enhancement of rj^ in contrast leads to slower relaxation 
which asymptotically follows a power law decay t^". For this 
reason, to secure the occurrence of normal transport the condition 
Vjyj.'jj » Trei must be well assured, with time scales separated by 
several orders of magnitude. Since the potential curvature can be 
estimated as K = (2it / L)^ AV , and the anomalous relaxation rate as 
[53] Trer(f;c/K)'^", one can see that the ratio ly^/AK is important 
to determine the scale of Tjei. Even if the transport of a large cargo 
is anomalous for AF«0.47 eV in Figure 3, the reduction of rj^ 
rapidly enforces normal transport, similarly to the reduction of 
Vtum- Therefore, for realistic turnover frequencies the molecular 
transport by kinesins is expected to be normal for vesicles of a 
typical radius of 30 nm, and possibly up to 100 nm. Anomalous 
transport emerges for above cargos with radius -300 nm or larger. 
This explains why the same molecular motors can mediate both 
normal and anomalous transport in living cells depending on the 
cargo size. The occurrence of an anomalous transport exponent is 
thus reconciled with the normal transport behavior for small cargo 
at lower turnover frequencies. 

Our research provokes a number of follow-up questions. Thus, 
what happens if we relaxed the assumption of a rigid motor-cargo 
linker molecule? In that case, the large subdififusing cargo is 
elastically coupled to a molecular motor, that possibly still operates 
normally in the absence of cargo. We are currently investigating 
this generalization for realistic spring constants of the linker. 
However, qualitatively the results remain very similar. Another 
question is prompted by the experimental results in Ref [33], 
suggesting that the motors can collectively transport several 
magnetosomes jointed into a chain. A generalization of our non- 
Markovian model to such collectively operating motors would be 
important for our understanding of large-cargo transport in living 
cells. 

Our model presents a good starting point for future research 
and further generalizations. Understanding how molecular motors 
perform in the viscoelastic cytosol of living cells despite the 
subdiffusion of the free cargo is compelling. Our findings open 
new vistas to the old problem of intracellular trafficking, 
reconciling seemingly conflicting results for the motor-cargo 
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dynamics under different conditions. Finally, our results will be of 
crucial importance for the design of new technologies of motor- 
driven particles and drug delivery in the crowded cytosol of cells. 
We are confident that our findings will prompt a series of new 
experiments on the dynamics of molecular motors under realistic 
conditions in living cells. 

Methods 

The numerical approach to integrate the generalized Langevin 
equation (2) with a power-law memory kernel rests on the 
approximation of the memory kernel by a sum of exponentials 
[34,53], 

JV 

1=1 

The rate constants v, and elastic constants ki are chosen to obey 
a fractal scaling [54], v, = vo/6'^', A:,ocv", with a dilation 
parameter &> 1, and a shortest memory time v^' in the hierarchy. 
Due to the scaling property rj^{ht,V(i) = h^"r]^{t,hvQ) this choice 
indeed provides a power law regime, )7„(?)ccf^°' on time scales 
Vq"' « ?«Tn]ax = 6'^^V^'0j with small superimposed logarithmic 
oscillations. Physically, this corresponds to representing a visco- 
elastic environment by auxiliary Brownian quasi-particles with 
coordinates x,-. They are coupled to the central Brownian particle 
by elastic constants kj and subjected to the thermal noises and 
frictional forces with viscous frictional constants =fc,/v, [34], 

N 

noX=f{x,t)- ^fc,(x-Xi) + ^o(0 > 
i= 1 

t]iXi = kiix - Xi) + ^/2r]ikBTi,i{t) , (4) 

where f{x,t)= —8V(x,t)/dx—fo, and are uncorrelated 

white Gaussian noises of unit intensity, <[ij{t')^j{t)y = dijd{t — t'), 
which are also uncorrelated to <^o(0- To have a complete 
equivalence with equations (2), (3), the initial positions x,-(0) are 
sampled from a Gaussian distribution centered around x(0), 
<x,(0)> = x(0) with variance <[xi(0)-x(0)]^> = A:^r/A;,. We set 

fc- = vo'/efr^(,_i)„j^(i_„)_jj (5) 

and use b=lO which leads to a maximal relative error with respect 
to the exact power-law of less than 4% [50], for a = 0.5, on 
relevant intermediate time scale. The effective relative friction 
coefficient ileti = *lefl/lO'^^ used as a parameter in our 
simulations. The parameter Vq controls the matching accuracy of 
our model with the model of fractional Langevin dynamics (the 
memory kernel )7„,(?) = '?3? "/r(l — a)) at short to intermediate 

times, where initially for times f«Tin = (f7o/'/a)''^*' °'' [51] the free 
diffusion is normal, (^6x^(t)y~2Dot with DQ = kBT/r]Q in accor- 
dance with the Einstein-Stokes relation. At longer times, 
anomalous diffusion emerges, <^Sx^{t)y~2Dxt^/r(l+oi) with 
Dai=kBT/ri^. The number of auxiliary quasi-particles controls 
the maximal range of subdiffusive dynamics, after which diffusion 
becomes again normal, <5x^(<))~2Z>efff, for <»Tmax = ^'^~ V^O 



with Deif=kBT/{rjgi[ + tj(,). The fractional friction coefficient is 

a-I / n Cai(h) 6' " , _jvn_5,), . 

'/a = '/effVaxA. where r= _ ^ [1 ^ >] is a. 

numerical coefficient of the order of unity, rxl.Ol for N>6, 
a = 0.5, and b= 10, with Ca(b)x 1.3 [53]. An interesting observa- 
tion is that in terms of Tmax and f/e(f, Tm~Tinax/»7eff ' Therefore, 
the effective viscosity jj^ff defines the time range of subdiffusion, 
from Tmax/'/eff' to Tmax> independently of b, and A''! For 
example, for »jeff = 10' ™d a = 0.5 one expects that subdiffusion 
win extend over 6 decades in time. In the simulations, we scale 
length in units of i- and time in units of Tq =-i-^»j/(47i;^i7o)> where 
Uq was taken equal to {70 = 0.25 eV (or about 10 ksT^oom with 
fixed temperature). The time step was (5/ = 0.01 and /niax was 
varied from 10* to lO'. The ratio U\/ Uf, was fixed to 0.65. Four 
different values of Uq were used: L^o = 1, 3/4, 2/3, and 1/2. The 
largest one corresponds to the largest potential barrier AF 
separating two potential periods of about 28 fc^Troom wO.7 eV, 
2/3 to about 18.7 A:£T'joom*0.47 eV, and the smallest one to 
14 ^B7^oom~0.35 eV. The larger Uq the more efficiently the 
motor works, as the probability of thermally activated backsteps is 
exponentially suppressed with AV/ksT. This is necessary to 
provide an ideal power-stroke operation at small turnover 
frequencies. However, the energy released from the hydrolysis of 
one ATP molecule wiU not be sufficient to perform one cycle of 
operation in the case AFq = 0.7 eV because of the energy derived 
from the hydrolysis of one ATP molecule does not exceed 0.62 eV. 
In such a case, to drive a cycle of two potential flashes one needs 
an input energy of at least 0.7 eV (0.35 eV per one conformational 
change, about one half of potential barrier separating potential 
periods, see upper inset in Figure 2). 

The power exponent of anomalous diflusion was fixed to a = 0.5 
to interpolate between a = 0.4 for an intact cytoskeleton [33,46] 
and a = 0.56 [33] when the actin filaments are disrupted. For the 
radius a = 300 nm of a single magnetic endosome [33] one 
estimates ijq ~5.65T0^' N-s/m in water of viscosity 
£=10^^ Pa-s. This yields to*2.24TO^^ s. Furthermore, we use 
Vi=V2 in the simulations. Hence, v = Vi/2, and the frequency 
vi=5T0^'* corresponds to v!*1.12 kHz. This is about the 
maximal possible turnover frequency for molecular motors [28]. 
For an effective particle ten times smaller, To«2.24TO^* s and 
Vl=5T0^' would correspond to the same physical turnover 
frequency. For a magnetic endosome we estimate '7eff = 10'', 
Vo = 0.1 and take N = S, which yields Tniax*22.4 sec and 
Z)c «3.68-10^ '* m^/s"'^, of the same order of magnitude as in 
experiments for a cargo consisting of several such endosomes [33]. 
Alternatively, with a smaller fi^^ = 10^ for vq = 1 and N = 9, is 
by a factor of ten larger, Z)aft!3.68T0^'^ m-^/s"'^, at the same 
Tmax- One obtains the same ZJj, by simultaneously increasing j/^fj- to 
j/cff = 10'' and Tmax to Tmax =2240 sec (iV= 11). 
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